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Abstract. The one-body and two-body density matrices in coordinate space and their Fourier transforms 
in momentum space are studied for a nucleus (a nonrelativistic, self-bound finite system). Unlike the usual 
procedure, suitable for infinite or externally bound systems, they are determined as expectation values of 
appropriate intrinsic operators, dependent on the relative coordinates and momenta (Jacobi variables) and 
acting on intrinsic wavefunctions of nuclear states. Thus, translational invariance (TI) is respected. When 
handling such intrinsic quantities, we use an algebraic technique based upon the Cartesian representation, 
in which the coordinate and momentum operators are linear combinations of the creation and annihilation 
operators a'^ and d for oscillator quanta. Each of the relevant multiplicative operators can then be reduced 
to the form: one exponential of the set {d^} times other exponential of the set {d}. In the course of such a 
normal-ordering procedure we offer a fresh look at the appearance of "Tassie-Barker" factors, and point out 
other model- independent results. The intrinsic wavefunction of the nucleus in its ground state is constructed 
from a nontranslationally-invariant (nTI) one via existing projection techniques. As an illustration, the 
one-body and two-body momentum distributions (MDs) for the ''He nucleus are calculated with the Slater 
determinant of the harmonic-oscillator model as the trial, nTI wavefunction. We find that the TI introduces 
important effects in the MDs. 



PACS. 21.60.-n Nuclear structure models and methods 
Nuclear reaction models and methods 
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1 Introduction 

The last few years the interest in the study of nuclei 
from both experimental and theoretical point of view has 
shifted from the investigation of one-body quantities (e.g., 
the elastic form factor F{q) and the momentum distribu- 
tion r]{p)) towards the investigation of two-body quanti- 
ties, with the aim of revealing more direct information on 
the dynamical correlations between the nucleons (short- 
range (SRC) and tensor). The two-body quantities are 
connected to the two-body density matrix (2DM) in co- 
ordinate or momentum space as are the one-body quan- 
tities connected to the one-body density matrix (IDM). 
The 2DM, besides being interesting in itself, allows the 
calculation of the expectation value of any two-body op- 
erator In addition to the 2DM p'^' (ri, ri/ , r2'), 
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we will also consider the two-body momentum distribu- 
tion (TBMD) r]^'^\p,k), which is the Fourier transform 
of the p'^' (ri, ri' , r2') in the variables ri — r2 and 
Ti' — T'2' and is connected to the two-nucleon spectral func- 
tion S{p, fe; E) via integration with respect to the energy 
E. 

A prominent role towards the experimental investiga- 
tion of the 2DM and related quantities is played by the 
study of the electromagnetically induced 2-nucleon emis- 
sion (7,NN), (e,e'NN) which can be carried out with high 
accuracy in photon facilities (Elsa, MAMI) and electron 
accelerators with high energy, 100% duty-cycle beams [Jef- 
ferson Lab, MAMI]. Past, present and near future experi- 
ments provide these useful data Theoretical 
methods to analyze the mechanisms of these reactions 
and to calculate the relevant nuclear two-body proper- 
ties are under continuous development. In particular, for 
the case of finite nuclei the generalized momentum dis- 
tribution ri{p,Q) T"^ , the two-body momentum distri- 
bution ry[^] (p, k) and other two-body distributions have 
been studied for Z ~ N, ^— closed nuclei, as well as the 
two-body density matrix for the nuclei ^He (31, "'^^O and 
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^°Ca and the two-nucleon spectral function S{p, k; E) 
for the nucleus ^^O (see [2] and refs. therein). 

One of the theoretical issues still under discussion is 
the proper consideration of the requirement of transla- 
tional invariance and therefore the conservation of the to- 
tal momentum of the system. The wavefunctions which 
have been used in the independent-particle shell model 
and in theories which take also dynamical correlations into 
account (e.g. Brueckner-Hartree-Fock, Variational Monte 
Carlo) satisfy the Pauli principle but not the translational 
invariance. As a consequence, they contain spurious com- 
ponents which result from the motion of the Center of 
Mass (CM) in a non free state. Effects from these (also 
known as CM correlations) are found in the calculation of 
almost every observable and make impossible the extrac- 
tion of information for the intrinsic properties of nuclei 
directly from the experimental data. In addition, there is 
ambiguity in the proper definition of translationally invari- 
ant operators which correspond to the different physical 
quantities. 

Many efforts have been made towards the solution of 
the CM problem. In some of them the treatment of the 
CM motion is built right into the theory one is using |12l 
I13II14 |. In the majority of the efforts the restoration of 
translational invariance is attempted after the wavefunc- 
tions have been developed. One such approach consisted 
of adding intuitively the three extra degrees of freedom 
of the CM to the ZA internal coordinates. Work along 
these lines using configurations from the harmonic oscil- 
lator model (HOM) has been carried out by Tassie and 
Barker and others. Most of the other approaches use 
projection techniques to define suitable intrinsic wavefunc- 
tions with coordinates referred to the CM. The pioneering 
works along this direction were made in the 50s by Garten- 
haus and Schwartz (GS) ^| and Peierls and J. Yoccoz 
(FY) ^21, with violations of the Gallilei invariance (GI), 
followed by Ernst, Shakin and Thaler (EST) ^M^ . with 
their critique on the GS transformation, Vincent [201 , She- 
beko et al pntS^ and others in the 70s [HII^ . Projection 
techniques have been proposed by Schmid and Griimmer 
24 , 25 at the beginning of the 90s and rather recently by 
Schmid and collaborators using harmonic-oscillator 
and spherical Hartree-Fock configurations [571 ■ We should 
add that Mihaila and Heisenberg have worked out the 
problem of CM corrections by expanding them as many- 
body operators |^ . It seems that if the wavefunction is 
very nearly factorable into a center-of-mass and an intrin- 
sic component all the approaches to treat the CM problem 
are equivalent, provided that the translationally invariant 
operators are used. Also, all the approaches can be carried 
out rather simply for the independent particle shell model 
with harmonic oscillator potential. 

The consideration of CM effects with one or more of 
the above mentioned methods have mostly addressed the 
light nucleus "^He and single-particle quantities: the ki- 
netic energy, the single - particle energies, the one-body 
density matrix, the matter and charge density, the elastic 
form factor, the dynamic structure factor, the momentum 
distribution and occupation probabilities, the one-body 



spectral function, the single-particle overlap function etc. 
There are few calculations for other light-medium nuclei 
such as i^C, i^N, 160, ^°C&. The consideration of CM cor- 
relations in two-body quantities has been limited so far, 
to the best of our knowledge to the potential energy. 

In this paper, we will start with the evaluation of the 
one-body density matrix in coordinate space and related 
one-body momentum distribution (OBMD). As mentioned 
above, CM effects have been considered in the calculation 
of IDM and OBMD before. By starting with such evalu- 
ation we want to present our method and technicalities. 
We then proceed to evaluate the two-body density matrix 
and two-body momentum distribution. We will take into 
account the CM correlations by using the EST prescrip- 
tion or fixed-CM approximation to construct the intrin- 
sic wavefunction. We also use a specific prescription for 
defining the corresponding intrinsic operators. As men- 
tioned before, the EST method has been introduced in 
refs. [THltT^ . Subsequently, it has been used in the cal- 
culation of the elastic and dynamic form factor and the 
momentum distribution of ^He '22 , 29' . The EST intrinsic 
many-body wavefunction is constructed from a transla- 
tionally non-invariant one by projecting onto an eigen- 
state of total momentum using a non unitary operator 
which fixes the CM coordinate R to be equal to zero. This 
transformation has certain advantages compared to GS 
transformation |19j (in particular, it ensures a correct be- 
havior under Galilean transformation). It has turned out 
(see also "5?) that the GS transformation can be re- 
duced to the EST projection procedure. As for the relevant 
one- and two-body intrinsic operators, they are defined by 
replacing the coordinates and momenta by relative ones 
(Jacobi variables). Unlike the definitions of the overlap 
integrals with intrinsic wavefunctions used in ref. ji-iOj . 
we are dealing with the expectation values of the intrin- 
sic operators as they occur under the treatment of the 
aforementioned quantities. Subsequently, we present the 
way for calculating the matrix elements defined by the 
above intrinsic wavefunctions and operators using an al- 
gebraic technique introduced in ref. [23 and based upon 
the Cartesian representation of the coordinate and mo- 
mentum operators in terms of linear combinations of the 
creation and annihilation vector operators (d^ and a re- 
spectively) for oscillator quanta in the three different space 
directions. With this technique one avoids to deal with 
difhcult multiple integrals and therefore it seems to be 
the technique of choice in the case of systems with large 
number of bodies. The Cartesian representation is partic- 
ularly convenient in the case of wavefunctions constructed 
with Slater determinants. The application of the above to 
the evaluation of the OBMD leads to the derivation of 
the Tassie-Barker factor jj^l (exprgg^/4A, rp oscillator 
length parameter entering into the definition of d^ and 
d - see sect. 2) in a model-independent way. In addition, 
the evaluation of both OBMD and TBMD leads to other 
model- independent results. 

Next, the intrinsic OBMD and TBMD are evaluated 
in the independent particle shell model with harmonic os- 
cillator wavefunctions. Such a model leads to compact an- 
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alytical expressions. Moreover, it is expected that the re- 
sults derived with its use will be close to the ones that will 
be obtained with the use of more realistic single-particle 
wavefunctions (Woods-Saxon or Hartree-Fock), since the 
quantities under study are not defined in terms of the 
asymptotic behavior of the wavefunction for large r. The 
latter is wrong in the case of the wavefunction in the HOM. 

In this work, the above evaluation is carried out for 
the intrinsic OBMD and TBMD of the nucleus "He. It 
is expected that the CM effects are more pronounced for 
light nuclei. In addition, due to its high central density (al- 
most 3 times nuclear matter density) the nucleus ''He is 
a particularly appropriate system to search for the origin 
of SRC. The evaluation of the intrinsic OBMD of ''He has 
appeared before in the case of HOM "57 as well as other 
single-particle models 29 . The (not intrinsic) TBMD of 
"He has been studied in refs. |11II31| by including in the 
HOM Jastrow-type correlations via the lowest term of 
the so-called low-order approximation ,i32. , but ignoring 
CM correlations. A comparison of the present results with 
those of refs will reveal the relative importance of 

the CM and SRC correlations in the same nucleus. The 
effect of the different correlations is estimated by introduc- 
ing the quantity rj^'^^p, k)/r]{p)ri{k). We find that the CM 
correction reduces the width of OBMD and TBMD and 
in the case of the TBMD introduces a dependence on the 
angle between p and k, a shift of its peak in favor of op- 
posite momenta and significant deviations for large values 
of p and k and angles close to 180°. This last effect is also 
found when SRC are considered. Up to now, the TBMD 
of ^He has been experimentally studied at Jefferson Lab 



The paper is organized as follows. In Section 2, we 
present the general formalism of constructing appropri- 
ate wavefunctions that respect translational invariance us- 
ing the Ernst-Shakin-Thaler prescription and describe the 
evaluation of matrix elements of intrinsic operators usung 
the Cartesian representation. In Section 3 the definitions 
of the relevant operators for the intrinsic quantities under 
study, namely the one-body quantities IDM, form fac- 
tor and OBMD and the two-body quantities 2DM and 
TBMD, are introduced in terms of the relative coordi- 
nates and momenta (Jacobi variables). In Section 4 the 
above quantities are evaluated using the Cartesian rep- 
resentation. In Section 5 by considering specifically the 
independent-particle shell model with harmonic oscillator 
wavefunctions, results are derived and discussed for the 
OBMD and TBMD of the nucleus "He. Finally, in Section 
6 a summary of the results and hints for possible further 
work are given. 

2 Constructing intrinsic wavefunctions and 
matrix elements. The Cartesian 
representation 

Let us consider a nonrelativistic system composed of A 
particles (nucleons). The coordinate (momentum) vector 
of the a— th particle will be denoted by (pa)- Occasion- 



ally, we will use the generic symbol a, which may include 
spin and/or isospin degrees of freedom, but in most cases 
we will suppress these degrees of freedom for the sake of 
simplicity. 

In principle, the eigenvectors of the total Hamiltonian 
H of the system |lf>), which belong to the eigenvalue P 
of the total momentum operator P, can be written as the 
product 

|fp) = \P) Ifint). (1) 

Following ref. 18 , the bracket | ) is used to represent a 
vector in the space of the center-of-mass coordinates, so 
that P\P) = P\P). A ket (bra) with an index | • • •)« («(• • • |) 
will refer to the state of the a— th particle. The intrinsic 
wavefunction Wi^t depends upon the A—1 independent in- 
trinsic variables. These may be expressed in terms of the 
Jacobi coordinates 

1 " 

l„ = r-„+i - - ^ (a = l,2,...,^-l) (2) 

13=1 

or the corresponding canonically conjugate momenta 

1 " 
Va = ——riaPa+i-^Pp) (a = l,2,...,A-l). 

^ 13=1 

(3) 

The wavefunction ^p{ri, r2, ta) in the coordinate rep- 
resentation satisfies the requirement of translational in- 
variance, 

ifp(ri+a,r2+a, . . . ,rA+a) ^ exp{iP-a)'Fp{ri,r2, ■■■,rA), 

^4) 

for any arbitrary displacement a. 

When describing scattering processes, it is convenient 
to consider the initial target state |0) as a P-packet, 

|0) = J |<f'p>dP(ifp|0> = j c{P)\^p)d^P (5) 

(see also Ch. XI), with the normalization condition 

(0|0)=| |c(P)pd3p = l. (6) 

Being the exact iJ-eigenvectors, the states \i'p) belong 
simultaneously to the set of eigenvectors of the total mo- 
mentum operator P with eigenvalues P close to a given 
value Pt, e.g., Pt — 0. The final state of the recoil- 
ing nucleus is written in the form \^pi). Evidently, the 
wavepacket |0) is not translationally invariant. However, 
this shortcoming can be corrected by letting the width of 
the packet A go to zero at the end of the calculations, i.e., 
assuming that 

Hml \c{P)\^9{P)d?P^ 1 5{P - Pt)g{P)dP ^ g{Pt) 

for an arbitrary function g{P). 

This prescription has a transparent physical meaning 
being adequate to many scattering situations. With its 
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aid one can express the corresponding cross sections in 
terms of intrinsic quantities. Let us consider, for instance, 
the elastic scattering of a particle from the nucleus. In the 
Plane Wave Impulse Approximation (PWIA) and neglect- 
ing the Fermi-motion effects, the cross section of interest 
can be represented in the form 



(8) 



a{e) = lim / X(P')d3p'|('fp'|cxp[iq-fA]|0)P 

where q is the momentum transfer, 9 the scattering an- 
gle and K{P') the corresponding kinematical factor. By 
substituting |0) from eq. (0), using eq. (1) and writing 
= {^A — R) + R, we evaluate (tf'p'l exp[iq • r^JjO) = 
/ d^Pc{P){P'\ exp(iq.^)|P)(tf'i„t| exp [iq ■ (f ^ - RWint) 
J (fPc{P)SiP' P q)Fint(q), that leads to 



a{9)^K{Pt + q) Fint(g) 



where the elastic (intrinsic) form factor (FF) Fint(q) is de- 
termined by 

Fintiq) = (tf'int|exp[iq-(fA-^)]|lp'int) = (!fi„t|Fint(9)|^'int)- 

(9) 

Use has been made, on the one hand, of the fact that in- 
trinsic and CM operators commute with each other and, 
on the other hand, of the condition ((T)). The Fint{q) is typ- 
ical of the quantities of interest, viz., it is the expectation 
value of an operator that depends on intrinsic coordinates, 
namely — R- 

As we have mentioned in the Introduction, our aim is 
to calculate expectation values of one- and two-body oper- 
ators in the nuclear ground state (g.s.) taking into account 
the requirement for TI. We are interested, in particular, in 
intrinsic quantities which appear in analytical expressions 
describing various scattering cross sections off a nucleus 
(in general, a finite system). Such quantities include, be- 
sides the elastic FF F{q), the particle density p(r), the 
dynamical FF S{q,oj), the OBMD ri{p) which is often as- 
sociated with the one-body spectral function P{p, E), and 
the TBMD ry'^l (p, fc) that is related to the two-body spec- 
tral function S{p,k;E). 

At the initial stage of the calculations, the nuclear g.s. 
is represented by a P - packet as introduced above. Then 
the main task is to construct the TI wavefunctions jtf'p) 
in a tractable manner, so that the CM-motion separation 
can be achieved. It is important to properly define the 
quantities of interest in terms of intrinsic coordinates, as 
was already done for the FF in eq. iPJ. We will tackle this 
issue in the next section. 

First of all, let us consider a Slater determinant 

|Dct> = ^ E ep^{l0pi(l))l</'p.(2))---|0p.(A))} 



VESa 

(10) 

as the total wavefunction |tfo) for an approximate and con- 
venient description of the nuclear g.s., in the framework 
of the independent-particle model or the Hartree-Fock ap- 
proach. In eq. (|10() . e-p is the parity factor for the permu- 
tation V, (pp^ are the occupied single-particle orbitals and 



the summation runs over all permutations of the symmet- 
ric group Sa- 

The wavefunction H10(l exemplifies wavefunctions which 
do not possess the property of TI, eq. Q). Obviously, any 
wavefunction that is constructed by acting on |Det) with a 
two- or a three- body correlation operator (e.g. a Jastrow 
correlation factor) will not be translationally invariant ei- 
ther. There are different ways to restore TI if one starts 
with a "bad" wavefunction such as |Det) [T7lll8n25l 
I26| . We shall employ the so-called "fixed-CM approxima- 
tion", or EST prescription However, other projection 
recipes can be applied without essential changes - we will 
come back to this point at the end of the Sec. 4. Within the 
EST approach, the approximate complete wavefunction is 
determined by 



(11) 



and the intrinsic unit-normalized wavefunction is given by 



EST\ 



{R^m/mmmY^^ (12) 



where {R — 0| is the eigenvector of the CM operator, R — 
J2a=i ^ct- have used the relation \R — X){R — 
X\ = S{R — X). For simplicity, we confine ourselves to 
the case of identical particles with mass m. The complete 
EST wavefunction can now be represented as 



1^W^)=Upm/m{2nfS{R)\'F) 



1/2 



(13) 



where, following |18| . we have introduced the projection 
operator {Up = Up) 



L/p = (2^)3/2|P)(fi = 0|. 



(14) 



Given an operator A, its matrix elements with the TI sym- 
metry can be written in the form 



{^'\u+Aup\^) 



[(i^"|(2^)3^(P)|if")('^l(27r)35(i2)|tf')]i/2 

(15) 

Its expectation value in the g.s. 



|0> 



d'Pc{P)Up\%)/[{% I i27TfSiR)\%)]'/^ 



is expressed in terms of the expectation value of the oper- 
ator U^,AUp, 

f f ^ , , {%\u+Aup\Wo) 

{o\A\o) ^ d'p dW(POc(p) ; / ^ 

J J (tf'o|(27r)3(5(P)|ip'o) 

(16) 

In addition, if A is an intrinsic operator Ai^t, acting only 
on the space of intrinsic variables, we find 

AeST = (0|iint|0) = {%\AEST\%)/{%\SiR)\%) (17) 
AeST =\R= Q)Ant{R = 0| = 6{R)Ant - iint<5(H). 

(18) 
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When deriving eq. (|17l ) we have employed the relation 
{P\P') = S{P - P') and equation ® . 

It has been shown 22 that the calculation of expecta- 
tion values of many-body operators like Aest can be sub- 
stantially simplified using the Cartesian representation. 
In this representation the coordinate (momentum) opera- 
tor fa (p^) of the a-th particle is expressed through the 
Cartesian creation and annihilation operators and d , 

f = -^(a+ + a) p = i^(a+-a) r^po = I, (19) 



obeying the commutation relations 
, [aua^. 



(20) 



which are the stepping stones in what follows. The indices 
j = 1,2,3 label the three Cartesian axes x, y, z. 
As the "length" parameter tq one can choose the os- 
cillator parameter of a suitable harmonic oscillator basis 
in which the nuclear wavefunction is expanded. Its basis 
vectors InxnyUz)!® . . ■'Si\nxnynz)A, where the quantum 
numbers rij,, Uy, take the values 0, 1, ... , are composed 
of the single-particle states 

\n,nyn,) = [njuyln^fi [a+]"^ [a+j"" [a+] |0 0), 

(21) 

which are the eigenstates of the Hamiltonian iJosc = w(d^ 

2 ' 



a+l) 



Hosc\nx riy Uz) = {ux + Uy + + |) w \nx Uy n^) , 

where uj is the oscillation frequency along the three axes 
a;, y and z. We use the system of units with h = c = 1. The 
single-particle wavefunction in coordinate representation 
is written 

(r I Uy Uz) = tpn^ ix)lpny iy)'<Pn, (z) , 

where |SIES1 

V'„(s) = [\/^2"n!ro]"*iI„(s/ro)exp(-sV2r2) 

and Hn{x) is a Hermite polynomial. By definition, the 
oscillator parameter equals tq = [mu] ~ 2 . 

The general idea in subsequent manipulations is to 
bring a given operator into a form with normal order- 
ing, in which the destruction operators d are to the right 
with respect to the creation operators (see Sees. 4, 5). 
For this purpose, we will also make use of the operator 
identity 

(22) 



which is valid for arbitrary operators A and B if the opera- 
tor C ~ [A, iJ] commutes with each of them. In particular. 



3 The intrinsic density matrices and related 
quantities 

In the preceding discussion it is implied that the operators 
of interest have been expressed in terms of the relevant 
coordinates, e.g., intrinsic ones. It is not always straight- 
forward how to do this. Here we refer mainly to the defini- 
tions of n— body density matrices (nDM's). For instance, 
it is a common practice |36II37II38| . to write the IDM in 
coordinate representation as the expectation value 



pW(r,r') = A(if|/5[il(r,r')|tf') 



(24) 



of the projection operator /5[^l(r,r'), 
pW(r,r') = \v)AA{r'\ 

Iri . . . rA-\r)Ari_ . . . drA-i(ri . . . rA-\r'\ 



\ (23) 



= exp(-ip^ • r)\rA = 0){rA = 0| exp(ip^ • r') 
= cxp(-ip^ • r)5{rA) exp(ip^ • r') 

in a given unit-normalized state ^ . Its diagonal elements 
give the one-body density distribution p{r) — p^^\r^r). 
The off-diagonal elements p^^\r, r') provide a measure of 
the correlation between the probabilities to find a particle 
in the two positions r and r' while all the other particles 
are kept fixed. Such a definition seems to be satisfactory 
in the case of infinite systems, or systems bound by an 
external potential, e.g. the electrons of an atom. However, 
it is apparently problematic for finite self-bound systems 
like nuclei, where the constituent particles are localized 
around their CM due to their interaction. Therefore, we 
prefer to deal with the intrinsic particle distributions that 
depend only on intrinsic wavefunctions and Jacobi coor- 
dinates. Only such quantities are of physical meaning in 
the case of finite self-bound nonrelativistic systems. In the 
next subsections this will be demonstrated for the intrinsic 
IDM and 2DM and related quantities. 



3.1 The intrinsic one-body density matrix and 
momentum distribution 

The intrinsic IDM in coordinate space may be defined as 

pW(r,r') ^A(tf'i„tlA™t(^,^')|S'int) (25) 
- = r)(^^_i - r'lS'i^t) (26) 

= A . . . d''U-2l'lMi, r) 

xtf'int(6,.. .,1^-2,0, (27) 

so that the normalization condition J d'^rp^^^^{r,r) = A is 
satisfied. We would like to emphasize that this is not an 
"imposed" definition. It appears naturally when evaluat- 
ing the dynamical FF [23 (or its diagonal part, if one uses 
the terminology adopted in Chapter XI of the monograph 
[331), which is related to the intrinsic OBMD ^ 



if [Al, Bj] — CSij for ~ 1, 2, 3) and C is a c— number. 



?7int(p) = A(^int\fHntip)\1'in 



(28) 
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with if the relations — P jA — r\^_^ and ~ Pa — 

mM = Sip - Pa + P/A) = Sip - f,,_,) (29) ^"'^^ conrpleteness of the 

^\Va-i=p){Va-i=p\- (30) 

[11 , „ i^intfAi, Aa) = / d^xd^yd^x'd^y'6(x + Ai - x') 

The OBMD is the Fourier transform of the IDMpj^'t (r,r'), J 

- (2.)- j dVdV exp [ip . (r - .')]^l(r, r'). ^ ^, ^3^^ 

(31) . 

See also ref. 40' . At the same time, the intrinsic one-body with 

density Pint{f) is the Fourier transform of the elastic FF ^p] , ' '^ it \/t 'i 

determined by eq. ©, or inversely, Pint(^' y;x ,y ) = \^a-2 = (£a-2 = ^ \ 

®\s.A-i^y){iA-i = y'\ ■ 

Fint{q)^jJ c"'-^PnAr)d\. (32) 



The latter is the intrinsic 2DM operator in coordinate 
space. It follows from eq. (|38ll that 

fj^{p,k) = (27r)-^y d'xd^yd''x'd'y'f^^{x,y;x',y'). 



From eq. ^ it follows that pi„t(r-) = A(ifint|Ant(r)|!Z'int), 
where 

An.(r-) = ^(r- - + ^) ^ Sir - A^L^,). (33) ^ ^^^^.^^ ^ , _ 

We notice that ^ g^p^i^, . _ y/)] (39) 

r n 3 

Pint(r') = -j^ P\^nti^~['^' ^T'^)- (34) Unlike the usual relationship 

In other words, the intrinsic IDM does not have the prop- (p, k) = fp'^ {p, k; p, k) 

erty p^^\r) = p[il(r,r) which can be justified for infinite _ /" 3 , , 3 , ip.(,_,') 

systems, although it has often been exploited in approxi- ~ v^^/ a ra sa r a s e e 

mate treatments of finite systems (cf., however, ref. -3(5, ^j2]/ / /\ 

where an alternative definition of the IDM for finite self- i*", , ), 

bound systems was proposed). where ryl^l (p, k) is the TBMD operator and pP^ (r, s; r' , s') 

{fj^^\p,k;p',k')) the 2DM operator in coordinate (mo- 

3.2 The intrinsic two-body density matrix and mentum) space as defined, for example, in refs. EHEIl, 

two-body momentum distribution the r.h.s. of eq. ^ contains a shift k/{A - 1) of the 

argument p, which may be negligibly small when the par- 
The formulation presented above will now be applied to tide number A increases. However, this is not the case for 
the 2DM. In particular, we will focus on the TBMD, usu- few-body systems, 
ally defined as the diagonal part of the 2DM in momen- 
tum space |11II31| . As we have already discussed, the rele- 
vant definitions require some revision in the case of finite, 4 The intrinsic density matrices and related 
self-bound systems. Here we will consider the expectation quantities in the Cartesian representation 
value 

[2], The intrinsic quantities are defined above in terms of op- 

^inUP- f') = ^(^ - ^{'i'intmPA-i -A^~P) erators, which can be written as products of A operators 

X(5(p^ — — fc)|tZ/int) acting on the subspaces of the separate A particles. Here 

ni we will show their evaluation in the Cartesian representa- 

^yl(A-l)(tf'i„tlCl(P,fc)l'^int), (35) tion. 

that can be interpreted as the TBMD with respect to the daltk^FF 'f' S""' ^^^"'^ ^^"^ operator related 

intrinsic momentum variables. We can write for the oper- ^ e e as ic , eq. 



[2], 

'"t'^^'"'^ Fi„t(9) =exp[ig-(fA-J«)] ^e"'"*-'' p-'^x^-'Jp'x±T'-A-q 



ator 7?i'„i(p,fc) ^ , , ... .-1 „ .-A-i 



41{P,k) = C^^y^ J d^Aid3A2e-'P-^ie-"=-^^i;i„t(Ai,A2). Notice that 

/5i„t(r) - S{rA -R-r) = {2^)'^ J e~"''^Fint{q)d\. 

Using eqs. (O, lEOl and we find 
-Eint(Ai,A2) = exppAi •T7^_2]exp[i(A2 - ^Ai) •r)^_i], , + . 

(37) e 'a 1 — Q s/2a'^ ^ ^ ' = e TA^e v^-* e v^-* 



(36) 

The operator i?int(Ai, A2) is expressed in terms of the Ja 
cobi variables, 
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-r-q 



n 2 2 

T2 4 , 



V2 



q-a ' 1- 

e 



q-a 



Oi(a)...OA-i(a)OA(/3) 



whence 

^'int(g) 

where 

" = -i7fe9 ' f3 = i^^il-A-')q. (40) 
Henceforth we use the notation 

0,(2) = e-^*-°^e^-"% (41) 



with i = 1, 2, . . . , A and any complex vector z, for a re- 
curring operator structure. In fact, one finds that all the 
intrinsic one- and two-particle operators of interest con- 
tain this operator structure: 

Oiiz) ■ ■ ■ dAMz)6A-lix2)dA{xi), 

where the vectors z,X2,Xi are related by the equation 
{A—2)z+X2+Xi = 0. In the case of one-particle operators, 

X2 = Z. 

Therefore, according to eq. ©, the elastic FF can be 
written as 



^int(q) 



■(<Z'int|Oi(a) . . . OA-iia)OAm^int) 

(42) 

where the vectors a and f3 are determined by eqs. 14UII . In 
the r.h.s. of the equation we find the Tassie-Barker (TB) 
factor exp{rQq^/AA) 15 with the "length" parameter tq. 
In our approach this factor results from the specific struc- 
ture of the intrinsic operator i^i,it('7), being independent 
of the nuclear structure (in general, the structure of the 
finite system under study). As it is well known, the TB 
factor appears directly in calculations, where the nuclear 
ground state is described by the simple harmonic oscillator 
model. 

In the fixed-CM approximation, as we demonstrated 
in Sec. 2 by means of eqs. (fTTji . (fH^ . one has to evaluate 
the ratio 

fkst(.) - (43) 



with 



FESTiq) = <5(H)Fint(q). 



Using the integral representation {2Tr)^6{R) = / exp(i^ • 
^)d'^^ and applying the same technique as before, one can 
show that 



i^EST(g) = c-(i--)— J d^fe" 

xdi{a')...dA-i{a')dAm 



(44) 



where 



'V2A 



/3'-i^K+(l-A-l)q] 



When calculating expectation values like those in eq. H43|) , 
the representation 144|) is especially helpful if the wave- 
function Wq is a Slater determinant or a linear combina- 
tion of Slater determinants. This has been demonstrated 
in particular in refs. |22ll41j . where the single-particle or- 
bitals entering the Slater determinant are eigenfunctions 
of a harmonic-oscillator potential with oscillator param- 
eter tq. Recently, similar calculations have been carried 
out beyond the HOM with the single-particle orbitals 
approximated by a truncated expansion in the Cartesian 
basis vectors of eq. (|21|l . 

The intrinsic IDM operator can also be expressed in 
terms of the Cartesian operators and a. We can rewrite 
the operator p\^tir,r') from eq. as 



A'int 



{r,r') = e-'^-*-i''5(|^„i)e'^^-i''' 
= (2^)"'/ d3Aexp[-i(p^- 
X exp[i^A • {r-A - R)] 
X exp[i(p^ - f ) • r'] . 



A) 



(45) 



As before, the general idea is to bring this operator in 
a form with normal ordering, where the destruction op- 
erators CLA are to the right with respect to the creation 
operators d\. To do this, we note that 



' exp 



• {r-A - R) 



exp[ 



X exp[i 



A-1 4 
A To 



A~l y2A 



^^Mi^,^^-{aA-^)], 

where £) = J2a=i is the "collective" destruction op- 
erator with the property [Di,Dj'] = ASij — 1,2,3). 
After some modest effort we arrive at the following result: 



Pint 



(r,r') = (27r)"^exp[- 



d^Aexph^^A^lexphif • (r + r')] 



X exp[(^(r - r') + i^.^X) ■ (a+ - ^)[ 
xexp[(-^(r-r')+iA:^ A) •("-§)] 



A-l V2A' 



The operator PestI'"'^'''') — ^(-^)p!nt(''' '''') '^^'^ be repre- 
sented in a similar way. Furthermore, taking into account 
eq. H31|) , the corresponding OBMD operator fjint [p] equals 



xexp[^y.(a+-^)] 
xexp[^y.(a-^)] (46) 

(2^)-3y d3yexp[zp-y]exp[-^4y'] 
x6i(x)...6a-i(x)Oa(7), (47) 
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As for the operator ?7est(p), we find x di{C') • . . Oa-2(C')Oa-i(72)Oa(7'i) , (53) 

/2 where 

X eM-V4e] exp[iy ■ p] I2 = + ^(^^ + ^A) 

xd,ix')...dA-i{x')OAil'), (48) and 

where 

j = -£^[l-^)y ^ ^ (49) The respective TBMD can be written as 

V = i;ii^-^(i-i)y , x'-i^i+^y-m 

With the help of eqs. lUZI), lUHll, (EHl, (EOl, the OBMD 



VESTiP^k) = (27r)-M(A-l) J d3yld^Ae-'P-(^+^/2) 

^^^-i..U^A/2) ^(^ + AAA-A/2) 

. ^ ^ (^ob)EST(p^ ^^"^'"^ 
"^""'^^^^ ('Z'o|<^(fi)|«'o) with 

N{\i,X2) - (<fo|^EST(Ai,A2)|<fo)- 

can be evaluated like the elastic FF i^EST(g) in the fixed- Qne can verify that this distribution meets the sequential 
CM approximation. relation 
By comparing the expressions 

r 1„ [rj^^{p,k)d'p^{A~l)rjEST{k). (55) 

'7int(p) = (27r) dy exp[iy • p] exp[iy • (p^ - ^4-^)]' 

The exposed method can be helpful in more general 
following from eq. 129|) , and situations when one has to handle translationally invariant 



^int(y) = exp[iy • (rA - R)], 



states of the kind 

\'I'^) = \P)W^nt), (56) 



we realize that it is sufficient to apply our algebraic tech- 
nique to the operator with the normalized state vector in the (3A— 3)— dimensional 

intrinsic Hilbert space \^^^) defined by 

exp[iy ■ {bi{a\- -ID^) + b2{aA- ^D))]- fGl^n) 

. „ ' [(•^o|G)(G|vI/o)]i/2- 
Then we can generate the intrinsic operators i^int (or pint) 

and 77int by changing the values of 61, 62. Here, |G) is any arbitrary vector in the CM space so that 

Similar manipulations lead to the following expression the scalar product (Gj^'o) represents integration of the 

for the operator Eint that enters in the definition of the CM variable only. Such general cases were considered in 

TBMD operator ??|'l(p,fc) (see eqs. (EBJ-lEHJ): ^ ^^"^ expectation value of the intrinsic operator 

iint in the state |0g) = / d3pc(P)|P) \^^^,) is 



2 \ 2 2.2 

Pn^ A-2 vj^ 



i;i„t(Ai, A2) =e 8 e 

x6i(C) • • • Oa-2(C)Oa-i(72)Oa(7i),(51) 

where 



Ag = (OcliintlOG) = (>Z'o|iG|!f'o)/(<f'o|G)(G|*o), (58) 



Ag - |G)Aint(G| - |G)(G|Aint = Aint|G)(G|. (59) 
(Cf. eqs. (24)-(25) in ref. [El.) 
71 = — ^A) 72 = ^={^-^A + iA) , The operator Aq, can be represented in the two equiv- 



alent forms 



' A ' /i__/j3D / ^3 jjii ry* r TD' 

We have set A — (Ai 4- A2)/2 and A = Ai — A2. In order 
to obtain the TBMD operator in the fixed-CM approxi- 
mation, one first needs to evaluate 



Ag = j d^P J d^P'G*{P')G{P)M{P,P')Aint (60) 



and 



£;est(Ai, A2) ^ (2nfS{R)E,^,i\,,X2). (52) = J d^R J d'R'G*{R')G{R)M{R,R')A,,, (61) 

Again, after some algebra one can show that ^^^j^ projection operators M(P, P') = |P)(P'| and 

. r ^ M{R,R') = |P)(P'|. The functions G(P) = (P|G) and 

^est(Ai,A2) = d ne " e ^0 ^ e t^^" q(^ji^ ^ (jj|q) represent the vector |G) in the two bases 
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composed of the P - eigenvectors (i.e., P\P) — P\P)) 
and R - eigenvectors (i.e., ii\R) = respectively. 

So far, the function G has been totally arbitrary. If 
we set G{R) = S{R), i.e., |G) = |H = 0) and G(P) = 
(27r)~^/^, we repeat the EST prescription with the intrin- 
sic state |!f'i^t''')j ^l- dJ' ^^"i tlis expectation value: 



iEST 



{%\S{R)Aint\% 

{%\S{R)\%) 



(62) 



If, on the other hand, we set G{P) = 5{P), i.e., |G) = 
|H = 0) and G{R) = (27r)~^/^, we arrive at the no-fixed- 
CM approximation after Peierls and Yoccoz J7] in con- 
structing the intrinsic state 



PY ^ jP = 0\%) 

[{%\SiP)\Wo)V/^' 



(63) 



The expectation value of interest is 



A 



PY 



{%\6iP)Aint\^o) 
i^^olSiPWo) 



(64) 



Note that within the PY approach the TI state vector for 
a nucleus moving with momentum P is approximated by 
jif'p"^) = \P){P\'Pq). Such a model is incompatible with 
the Gallileo invariance requirement. The restoration of 
Gallileo invariance was discussed by Peierls and Thouless 
j42| . The unit-normalized vector given by eq. (|63|l corre- 
sponds to the PY prescription for the nucleus rest frame 
only. In this context we would like to refer again to the sev- 
eral calculations in nuclei with the ansatz H63fl performed 
recently by Schmid and his colleagues |ESII77j . 

Previous experience prompts us to unify forthcoming 
calculations of quantities like (|62|l and H64|) by introducing 

one i5— function, S{ciD^ + C2-D), where ci, C2 c-number 
parameters, instead of the two functions S{R.) and S{P). 
Along these guidelines, one is able to build up a family 
of generating functions which can be calculated with the 
help of the algebraic technique developed in this article. 
Finally, the following representation 

M{R, R') = exp(-iP • R)6{R) exp(iP • R') 

might be useful in calculations employing a refined de- 
termination of the weight function G. In particular, as 
outlined in ref. JS] , the "best" G must be chosen so as to 
minimize the expectation value of the intrinsic Hamilto- 
nian. See also ref. [201 1 devoted to an optimal separation 
of CM motion in many-body systems. 

Closing this section, we should note that within the 
simple HOM, where the single-particle states are described 
as pure harmonic-oscillator wavefunctions, the total wave- 
function is always a product of the CM wavefunction and 
an intrinsic one. Therefore, for intrinsic quantities, one 
gets the same results with or without the use of a projec- 
tion technique. 



5 The intrinsic one-body and two-body 
momentum distributions: Application to "^He 
and discussion 

The many-particle operators encountered so far have much 
in common with each other owing to the operator struc- 
ture appearing in each of them, 

di{z) ■ ■ ■ dA-2iz)dA-l{x2)dA{xi). 

Let us derive their expectation value in the independent- 
particle model. The wavefunction is then a Slater deter- 
minant |Det) (eq. (|10|l ). The following formal result is 
straightforward: 



(Det|Oi(2) • • • dA^2{z)dA-i{x2)dA{xi)\Det) 
= (Det(-a;i, -a;2, -z)\'Det{xi,X2, z)), 



(65) 



where 



iDet 



X \c^,,_, {A - 2))e--'^-- {A - l))e--^- (A))}(66) 

The r.h.s. of eq. (|65|) can be represented as the sum of 
terms containing the matrix elements |Om l^^pj ) ihj,n^ — 
1,2, .. .^). Note that the permutations V are related to 
the single-particle states (not the nucleon labels). In the 
simplest case of the Os'* configuration, that we encounter 
in ^He nucleus, all we need is to evaluate the matrix 
element (0s|Om|0s). We have written \(j)os) = |Os)|(Tr), 
where |0s) and |(tt) are the space and spin-isospin parts, 
respectively, of \4>os)- In the HOM, (f>os{r) = (r|Os) = 
(-\/7r6)~^/^e~'' Z^** , where b — ro — l/po is the harmonic 
oscillator parameter. The state |0s) coincides with the 
lowest-energy state |000) which is the "vacuum" of the 
Cartesian representation, viz., d|000) = 0. As a conse- 
quence, the radial matrix element of interest is equal to 
unity and 



{0s^\Oiiz)O2{z)d3ix2)Oiixi)\0s^ 



1, 



(67) 



rendering the rest of the calculation trivial. We are now in 
a position to write down expressions for all the quantities 
considered, in the case of "^He. In particular, performing 
the necessary integrations and taking into account the nor- 
malization of each quantity, we get for the intrinsic OBMD 
and TBMD in the EST prescription 



'7est(p) = Vest{p) 



43/2453 



33/2^3/: 



-e 3 



p b 



[21 , 23/2 126^ 

%ST(P>fc) = — —3 — e 2^ 



(68) 
(69) 



Regarding the properties of these distributions, we find 
that they satisfy the relation H55|) : 



'yESTb' ^)<^'^P = 377EST(fc)- 



10 



A. Shebeko et al.: Density matrices of finite nuclei with an appropriate treatment of CM motion 



In the simple HOM (with the CM correction not included) 
we have EiniEI] 



463 

ViP) = V{P) = 3372*^ 



2 1^2 
~p 



TT' 

„2l.2 



(70) 
(71) 



We notice that, within the HOM, the expression for the 
OBMD in the EST approach, eq. (jSHl), is the same as the 
one without the CM fixation, eq. JTUJ, if one substitutes 
6 with •\/4/36. This results in a narrowing of distribution 
'/EST (p) /»7est (0) compared to r]{p)/r]{0). But this obser- 
vation does not mean that the CM corr ection reduces to 
the renormalization 6 bcM = a/ 4/36. In this connec- 
tion, let us recall the work [22], where the charge FF and 
the dynamic FF of ''He were calculated under some simpli- 
fying assumptions using the same oscillatory wavefunction 
and the EST prescription. Its main discovery was the nar- 
rowing of both the intrinsic density and OBMD due to the 
CM correction. In order to demonstrate this phenomenon, 
let us write the nucleon density in the EST approach. 



PEST{r) = 4- 



exp(-§T) 



^3/2^3 ' 



(72) 



that follows directly from previous derivations and rewrite 
eq. l(7n|) in the form : 



Po' 



(73) 



Here we have come back to the oscillator parameters rg 
and Po and introduced the new values fo = \ ^^^^''o and 



Po 



^ ^Pa- Meanwhile, we write i/^^ instead of 



3/4 to point out a trend in the A - dependence of this 
effect of center-of-mass motion. We notice that the corre- 
sponding quantities in the simple HOM 



p{r) 



r]{p) = 4 



TT^^'^Pa^ 



are obtained from eqs. H72|l - (|73|l when A oo . 

Thus we see (cf. j^) that the inclus ion of CM cor- 
rections gives rise to the two renormahzations, — s- fg = 

\J^^ro and po ^ po — \J^^Po, of the oscillator param- 
eter values, ro and po, in the density and momentum dis- 
tributions calculated within the simple HOM. Evidently, 
such changes are not accounted for by a hasty replace- 
ment of 6 by -\/4/36. (Therefore, the discussion in ref. 1J, 
p. 263, on the results of ref. j22j is incomplete.) One may 
say that we encounter a specific effect of shrinking of the 
density distribution p{r) and the momentum distribution 



7/(p). The term "shrinking" implies that each of these den- 
sities, after being CM corrected, increases in its central but 
decreases in its peripheral region. As it has been shown in 
the past |22II46| . such a simultaneous change of the one- 
body distributions plays an essential role in getting a fair 
treatment of the data on elastic and inelastic electron scat- 
tering off ''He. Notice that the product fopo = 1 — 1/A^ 1, 
unlike the relation ropo = 1. The commutation rules for 
the intrinsic coordinates r' = r R and conjugate mo- 
menta p' = p — P/A 

[ir% ip'),] = i%(l - 1/A) ; 1,2,3 

and the corresponding uncertainty principle are related 
to this deviation from unity |41| . Thus, the uncertainty 
principle is not contradicted by the simultaneous shrinking 
of the density and momentum distribution (see also 021, 
Lect. I, Suppl. C). 

A similar shrinking is encountered also in the two- 
dimensional surface given by the function 77^ggrp(p, fc) = 

^est(P' ^) O"^^ given by 777 '(p, k) = rp'^p, fc), at 

each value of the correlation angle 7. This effect is clearly 
visible if the TBMD of ''He is represented as 



A-1 



3/2 



77]^! (p, A:)e 



Now, in order to compare numerically calculations based 
upon the formulae (|69|) and (|71|l . we adopt, as in ref. [22] . 
the following way of determining the harmonic-oscillator 
parameter 6. We start with a general expression for the 
charge FF Fc{q) = fp{q)F{q) of a nucleus, where fp{q) the 
correction factor for the finite proton size. The factor fp{q) 
is approximated by a Gaussian, fp{q) = exp{—bp^q'^) and 
for the parameter bp we choose the value bp = 0.8 fm j43j . 
The FF of the ^He nucleus equals F{q) = exp(— ro^g^) in 
the simple HOM and F^sriq) = expl-f'^q'^) in the HOM 
with CM corrections. In each of the two models the value 
of the adjustable parameter tq — b will be determined 
by the requirement to reproduce the experimental value 
of the charge mean square radius of "^He, rims = 1-67 fm 
|77| . Then in the simple HOM (6 = 60) we have 



and therefore 60 = 1.197 fm. When CM motion corrections 
are taken into account (6 — 6cm), 



3A-1, 



''rms — o A ^CM 



2 A 

and therefore 6cm = ^4i6o = 1-382 fm. By adopting 
this fitting we therefore get the identical g-dependence, 
FEsriq) = F{q) = exp(— 60^9^), while the difference be- 
tween the respective OBMDs becomes more considerable. 
In fact, we have 



msTip) = 4 



A-l 



3 60' 



r3/2 



exp(- 



A~l 



pW) 
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4 

2 

1 

0.8 
0.6 
0.4 
0.2 


0.06 
0.05 
0.04 
0.03 
0.02 
0.01 




1 1 1 


1 1 1 

HO with CM corr. 




H01 




~\ H02 




\ P=0 


1 — ^-'"\"' 







1 1 1 


1 1 1 1 

HO with CM corr. 






; H01 






H02 






\ p=1 fm"^ 






1 r-O^r — ...J 



states) . Even for the finite non-interacting fermion system, 
^ = 1 — l/iyifp=:fc, (because rp^^{p,h) = ri{p)ri{k) — 
^\rii{p, k)\'^ holds). Deviations of £,{p,p) from this value 
show the effect of other-than-statistical correlations. For 
p k, deviations from 1 — 1/A is a measure of statistical 
and (or) dynamical correlations in a system of finite size. 
In addition, in the case of self-bound systems, deviations 
from this value account for the correlation due to the fixed 
center of mass of the system. 

The system of "^He in the simple HOM is a special case 
for which v = A and therefore 



1 1 1 


1 1 1 

HO with CM corr. 




HOI ~ 




H02 _ 




"--./^^ p=1.5fm"^ 




^-~4_ i- --r4-- 



-1.5 -1 -0.5 0.5 1 1.5 2 

Fig. 1. The TBMD of "^He for p || fc and p = 0, 1, 1.5 fm'S 
as a function of kp, where k = kpp, in the HOM model. Full 
line: with CM motion effects taken into account within the 
EST approach, eq. (I69L with 6 — 6cm; long-dashed line: in the 
simple HOM, eq. I)71|l . with b — bo (HOI); short-dashed line: in 
the simple HOM, eq. JTTJ, with b = bcM (H02), from refs. [TTl 



vs ^ 

The former differs from the latter by a substitution &o ~^ 

In Fig. 1 the TBMD of ''He is shown for p = 0, 1, 1.5 im-^ 
as a function of kp, where k is parallel to p and k — kpp, 
i.e., kp is positive (negative) for k in the same (opposite) 
direction as p. The TBMD that we have calculated us- 
ing the EST method, i.e., with CM motion corrections 
(6 = bcm), is plotted with full fines, while the TBMD 
within the simple HOM with h — ho \s plotted with long- 
dashed lines ("HOI") and the TBMD within the simple 
HOM with b — bcM is plotted with short-dashed lines 
("H02") taken from refs. HHEIl- One can observe the 
shift of the peak from kp = towards negative /cp's, for 
p ^ 0, due to the correlation induced by the fixed center 
of mass. 

In refs. |11II31| the dimensionless quantity 



^{p,k) EE ri^^\p,k)li^{p)'q{k) 



(74) 



was introduced, as a measure of correlations of statistical 
and dynamical origin as well as of finite-size effects. In 
the complete absence of correlations, ^ should be equal 
to 1 — 1/A. In the case of the infinitely extended ideal 
Fermi gas ^ is defined for p, fc < pp, where pf is the Fermi 
momentum, and if p 7^ fc, ^ = 1 (note that A — > 00), while 
for p = fc, ^ = 1 — (i' is the degeneracy of the particle 



i = l-l/v^ 0.75 



(75) 



for all p and k. The same does not hold after fixing the 
center of mass. Then we have: 



^est(p, fc) 



'yjs^SxlP' fe)/'7EST(p)?7EST(fe) 



= 0.89493e-sP"''"e"^'^"''"e-P'=''" (76) 

In Fig. [3 we plot ^ as a function of kp, where k — kpp, 
for p = 0,1,1.5 fm~^. In Fig. |3 log^o C is plotted as a 
function of cos 7, where 7 is the angle between p and k. 
The full, long-dashed and short-dashed lines correspond 
to the Cest(p, fe) of eq. (|76|l for p = fc = 1 fm~^, for 
p = 1, fc = 4 fm^^ and for p = fc = 4 fm~^, respectively, 
and the dashed lines to ^{p, k) — 0.75, eq. (|75|l . The effect 
of CM correlations is important. The EST TBMD favors 
momenta of opposite directions, as compared to the prod- 
uct of the two OBDM, while in forward angles ^est is 
significantly reduced. In refs. [1111^] the effect of short- 
range correlations (SRC) was investigated by including in 
the simple HOM Jastrow-type correlations in the calcula- 
tion of TBMD of the ''He nucleus using the lowest-order 
approximation of ref. |32| . Significant deviations from the 
independent-particle picture were found for large values 
of p and k close to 7 = 180° and 0° . In Fig. 3 the corre- 
sponding quantity ^ for the case p = fc = 4 fm~^ is plotted 
in logarithmic scale for comparison. It is anticipated that 
within the EST approach additional corrections due to 
SRC will appear mainly at high values of p and/or k and 
that they will be larger when p and k are antiparallel. 



6 Summary and Conclusions 

The intrinsic one-body and two-body density matrices in 
coordinate space and corresponding Fourier transforms in 
momentum space have been studied for a nucleus (a non- 
relativistic system) that consists of A nucleons (particles). 
We have seen how these quantities of primary concern can 
be expressed through expectation values of the A— particle 

multiplicative operators pfj^^J. and pj^j. sandwiched between 
intrinsic nuclear states. Our consideration is translation- 
ally invariant since the operators depend on the relative 
coordinates and momenta (Jacobi variables). To avoid a 
cumbersome multiple integration, we have developed an 
algebraic technique based upon the Cartesian representa- 
tion, in which the Jacobi variables are the linear combi- 
nations of the creation and destruction operators for os- 
cillator quanta in the three different space directions. In 
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Fig. 2. The quantity ^ of "He for p || fc and p = 0, 1, 1.5 fm"\ 
as a function of kp, where k — kpp, in the HOM model. Full 
line: with CM motion effects taken into account within the 
EST approach, eq. l|7()^ for b — Bcm; dashed line: in the simple 
HOM, eq. refs. 
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Fig. 3. The quantity logj^g ^ of *He for the indicated values of p 
and fc as a function of cos 7, where 7 is the angle between p and 
fc, in the HOM model. Full, long-dashed, short-dashed lines: 
with CM motion effects taken into account within the EST 
approach, eq. 1761 for b — bcM', dotted line: including short- 
range correlations within the LOA |11|. without CM motion 
effects taken into account; dot-dashed line: in the simple HOM, 
eq. 1751 (result independent of p, k and 7). 



the framework of the subsequent operations the normal 
ordering of the operators involved in pj^J. and pjj^jj. plays 
a central role in getting both the general results and the 
working formulae. 

The Cartesian representation is convenient and allows 
us to find simple links between the relevant distributions 
via the generating functions constructed here. In partic- 
ular, the OBMD rj{p) and the elastic FF F{q) can be 
deduced from one and the same generating function by 
changing the values of its arguments. 

In the course of such a procedure the so-called Tassie- 
Barker factors stem directly from the intrinsic operators 
(not the WFs). One should emphasize that these factors 
(different for different distributions) occur here reflecting 
the translationally invariant structure of the correspond- 
ing intrinsic operators. Each of them is a Gaussian whose 
behavior in the space of variables is governed by the size 
parameter ro (or its reciprocal po) and the particle num- 



ber A for a given system, but it docs not depend upon the 
choice of the g.s. WF. The latter can be a simple Slater 
determinant, include SRC or not, be CM-corrected or not, 
etc. In practical calculations, such WF's (in particular, the 
reference WF jtf') in eq. (|12|) ') are often expanded in the 
convenient HO basis functions. Therefore, in order to ex- 
ploit all the power of the HO algebra when manipulating 
the intrinsic operators, it is pertinent to set the working 
ro - value equal to the respective "optimum" value of the 
oscillator parameter. 

In order to realize all the general results, obtained 
above for the intrinsic density matrices and related distri- 
butions, the intrinsic wavefunctions of the nuclear ground 
state have been constructed using the prescription by Ernst, 
Shakin and Thaler, that leads to the so-called fixed-CM 
approximation. In this connection, we have demonstrated 
how one can unify the different approximate recipes of 
restoring the TI, if one starts with one of them, e.g., with 
the EST projection operator. As a specific example, ana- 
lytic expressions for the intrinsic OBMD r/(p) and TBMD 
7y[^l(p, fc) of the "'He nucleus have been derived within 
the context of the independent particle shell model, using 
harmonic-oscillator wavefunctions. When CM corrections 
are taken into account, the OBMD and TBMD are simul- 
taneously shrunk with respect to the nontranslationally 
invariant counterparts. In addition, the CM correlation 
introduces in the case of the TBMD a dependence on the 
angle between p and fc. A shift of its peak for p 7^ in 
favor of opposite momenta and significant deviations for 
large values of p and k and angles close to 180° are ob- 
served. The above calculation is relevant to the current 
experimental study of two-nucleon knock out off He iso- 
topes. For instance, we see similar behavior as the one 
found in the recent experimental study of the TBMD in 
^He 0. 

Of course, when increasing the momenta transferred 
to a residual system (in particular, under kinematic condi- 
tions where they are getting comparable with the nucleon 
mass) the corrections to TI breaking should be consid- 
ered along with incorporating relativistic effects such as 
the Lorentz contraction of nuclear WF's and a specific 
velocity-dependence of nuclear forces. The latter is typi- 
cal of relativistic one-particle-exchange models (see survey 
P5| and refs. therein), where, for instance, the nucleon- 
nucleon quasipotential, including recoil effects, is essen- 
tially nonlocal and prevents, in contrast to the nonrela- 
tivistic case, the separation of the CM motion of the rela- 
tivistic system (nucleus) from its internal motion. In other 
words, the corresponding four-momentum eigenstates can- 
not be factorized as a product of independent CM and 
intrinsic components. In this context, note a very instruc- 
tive work 0^ . Taking into account this distinctive feature 
of relativistic quantum mechanics, it is difficult a priori 
to say to what extent the approach developed here could 
be helpful in a covariant description of composite systems. 
Nevertheless, an encouraging example is found in Ref. [KT)] . 
where the fixed-CM approximation combined with an ap- 
propriate Lorentz contraction has been used for calcula- 
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tions of the nucleon electromagnetic FF's in the cloudy 
bag model with CM and recoil corrections. 

Using the techniques presented in this paper one could 
calculate the intrinsic TBMD of other than *He Z — N, 
£— closed nuclei within the context of the harmonic os- 
cillator model. In addition, within the present framework 
one could investigate the effects of short-range correla- 
tions on the intrinsic TBMD of "^He and other ^—closed 
nuclei, which are expected to be sizable for p,k > pp, 
by introducing Jastrow-type correlations. Also, other in- 
trinsic two-body quantities could be evaluated within the 
above general formalism (including other non-relativistic 
systems). Finally, the approach developed here may be 
helpful when evaluating the intrinsic one-body, two-body 
and more complicated density matrices in the HOM and 
in other independent-particle models (e.g., with single- 
particle wavefunctions of a potential well with finite depth, 
as it was shown in Ref. |41jV 
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